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In the past decade, the number of known binary near-Earth asteroids has more 
than quadrupled and the number of known large main belt asteroids with satellites 
has doubled. Half a dozen triple asteroids have been discovered, and the previously 
unrecognized populations of asteroid pairs and small main belt binaries have 
been identified. The current observational evidence confirms that small (<20 km) 
binaries form by rotational fission and establishes that the YORP effect powers the 
spin-up process. A unifying paradigm based on rotational fission and post-fission 
dynamics can explain the formation of small binaries, triples, and pair's. Large (>20 
km) binaries with small satellites are most likely created during large collisions. 


1. INTRODUCTION 
1.1. Motivation 

Multiple-asteroid systems are important be¬ 
cause they represent a sizable fraction of the aster¬ 
oid population and because they enable investiga¬ 
tions of a number of properties and processes that 
are often difficult to probe by other means. The 
binaries, triples, and pairs inform us about a great 
variety of asteroid attributes, including physical 
properties, composition, interior structure, for¬ 
mation processes, and evolutionary processes. 

Observations of binaries and triples provide the 
most powerful way of deriving reliable masses 
and densities for a large number of objects. The 
density measurements help us understand the 


composition and internal structure of minor plan¬ 
ets. Binary systems offer opportunities to mea¬ 
sure thermal and mechanical properties, which 
are generally poorly known. 

The binary and triple systems within near- 
Earth asteroids (NEAs), main belt asteroids 
(MBAs), and trans-Neptunian objects (TNOs) ex¬ 
hibit a variety of formation mechanisms (Merline 
et al. 2002c; Noll et al. 2008). As such, they pro¬ 
vide an invaluable window on accretional, colli- 
sional, tidal, and radiative processes that are criti¬ 
cal in planet formation. The distribution and con¬ 
figurations of the multiple-asteroid systems also 
provide a rich array of constraints on their envi¬ 
ronment, their formation, and their evolutionary 
pathways. 

Observations rely primarily on ground-based 
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telescopes and the Hubble Space Telescope (HST). 
For an up-to-date list of binaries and triples in the 
solar system, see Johnston (2014). We describe 
observational techniques only briefly because this 
material is available elsewhere (e.g., Merline et al. 
2002c). A few emerging techniques will be de¬ 
scribed in more detail. Likewise, we refer the 
reader to other texts for an extensive history of 
the field (e.g., Merline et al. 2002c) and highlight 
only a few of the developments here. 

1.2. History 

Early search programs for asteroid satellites 
were unsuccessful, returning negative or dubious 
results, such that the authors of the Asteroids II 
review chapter chose the prudent title “Do as¬ 
teroids have satellites?” (Weidenschilling et al. 
1989). The chapter provides an excellent discus¬ 
sion of the physics of several formation mecha¬ 
nisms that were postulated at the time. The per¬ 
spective changed with the flyby of (243) Ida by 
the Galileo spacecraft in 1993 and the discov¬ 
ery of its small satellite Dactyl (Chapman et al. 
1995; Belton et al. 1995). Ground-based efforts 
intensified and resulted in the discovery of a satel¬ 
lite around (45) Eugenia by Merline et al. (1999). 
Several other discoveries followed in rapid suc¬ 
cession. The relatively small sizes of the MBA 
satellites suggested formation in sub-catastrophic 
or catastrophic collisions (Durda 1996; Dores- 
soundiram et al. 1997). 

The discovery of MBA satellites, coupled with 
analysis of terrestrial doublet craters (Bottke and 
Melosh 1996a,b) and anomalous lightcurve ob¬ 
servations (Pravec and Hahn 1997), suggested the 
existence of binary asteroids in the near-Earth 
population as well. The unambiguous detection 
of five NEA binaries by radar cemented this find¬ 
ing and indicated that NEA satellites form by 
a spin-up and rotational fission process (Margot 
et al. 2002). Lightcurve observers reached the 
same conclusion independently (Pravec and Har¬ 
ris 2007). Both radar and lightcurve observations 
revealed that, far from being rare, binary asteroids 
are relatively common (Pravec et al. 1999; Margot 
et al. 2002; Pravec et al. 2006). By the time the 
Asteroids III review chapter was written, a more 


decisive title (“Asteroids do have satellites”) had 
become appropriate (Merline et al. 2002c). This 
review focuses on the developments that followed 
the publication of Asteroids III. 

1.3. Terminology 

Two- and three-component asteroids that are 
gravitationally bound will be referred to as bi¬ 
nary asteroids (or binaries) and triple asteroids 
(or triples), respectively. ( Triple is favored over 
the more directly analogous terms trinary and 
ternary because of long-established usage in as¬ 
tronomy). Asteroid pairs denote asteroid compo¬ 
nents that are genetically related but not gravita¬ 
tionally bound. Paired binaries or paired triples 
are asteroid pairs where the larger asteroid is itself 
a binary or triple asteroid. The larger component 
in binaries, triples, and pairs is referred to as the 
primary component or primary. The smaller com¬ 
ponent in binaries is referred to as the secondary 
component or secondary. 

There has been some confusion in the literature 
about the meaning of the word “asynchronous.” 
Here, we adopt the terminology proposed by 
Margot (2010) and later implemented by Jacob¬ 
son and Scheeres (2011b) and Fang and Margot 
(2012c). Binaries with an absence of spin-orbit 
synchronism are called asynchronous binaries. 
Binaries with a secondary spin period synchro¬ 
nized to the mutual orbit period are called syn¬ 
chronous binaries. Binaries with both primary 
and secondary spin periods synchronized to the 
mutual orbit period are called doubly synchronous 
binaries. If generalization to systems with more 
than one satellite is needed, we affix the terms 
synchronous and asynchronous to the satellites 
being considered. 

It is useful to present results for small and large 
asteroids. We place an approximate dividing line 
at the size at which objects are substantially af¬ 
fected by the YORP effect during their lifetime. 
For typical NEAs and MBAs, this dividing line 
corresponds to a diameter of about 20 km (Jacob¬ 
son et al. 2014a). We define very small asteroids 
as those with diameters of less than 200 m. This 
is the approximate size below which many aster¬ 
oids are observed to spin faster than the disruption 
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rate of a body with no shear or tensile strength 
ojd = where p is the density and G is 

the gravitational constant. 

We use two additional acronyms. The YORP 
effect is a radiation-powered rotational acceler¬ 
ation mechanism for small asteroids (Rubincam 
2000). The binary YORP (BYORP) effect is 
a radiation-powered acceleration mechanism that 
may expand or contract the orbits of some syn¬ 
chronous asteroids (Cuk and Burns 2005). 

2. OBSERVATIONS 

Several observational techniques are available 
for discovering, detecting, and studying binaries, 
triples, and pairs, each with its strengths and 
weaknesses. This section describes recent results 
and illustrates the complementarity of the obser¬ 
vational techniques that characterize individual 
asteroid systems and entire populations. 

2.1 Radar Observations of NEA Systems 

Radar has proven to be a powerful method 
of detecting secondaries to NEAs, enabling the 
discovery (as of September 2014) of the satel¬ 
lites in 71% of the 49 known multiple-component 
NEA systems, including 33 of 47 binaries and 
both undisputed triple systems. Of the 14 bi¬ 
nary NEAs discovered via optical lightcurve tech¬ 
niques, 6 have been confirmed with follow-up 
radar observations during later apparitions. Over¬ 
all, radar detections suggest that about one in six 
NEAs larger than 200 m in diameter are multiple- 
asteroid systems (Margot et al. 2002; Taylor et al. 
2012a), though 200 m is not a sharp cutoff. Three 
binary NEA systems identified by radar have pri¬ 
mary components with suggested diameters of 
120 m to 180 m: 2003 SS 84 (Nolan et al. 2003), 
(363599) 2004 FGn (Taylor et al. 2012c), and 
1994 CJi (Taylor et al. 2014). For comparison, 
the largest primaries of binary NEAs imaged with 
radar: (5143) Heracles (Taylor et al. 2012b), the 
possible triple (276049) 2002 CE 26 (Shepard et al. 
2006), and (285263) 1998 QE 2 (Springmann et al. 
2014), are more than an order of magnitude larger 
at >3 km in diameter. It is likely that ~8 km 
diameter (1866) Sisyphus has a secondary based 


on analysis of frequency-only observations ob¬ 
tained on four separate dates in 1985 (Ostro, pers. 
comm., 2001). 

Radar observations can be used to detect aster¬ 
oid satellites because of the ability to resolve the 
components of the system both spatially (along 
the observer’s line of sight) and in terms of fre¬ 
quency (due to Doppler shifts from the rotational 
and orbital line-of-sight velocities), resulting in 
a measurable separation between the components 
in two dimensions. Direct detection of a satel¬ 
lite in frequency-only spectra or radar images typ¬ 
ically occurs within one observing session and 
often within minutes of observation. The band¬ 
width of the echo of a component scales directly 
with the diameter and rotation rate. Thus, in 
a frequency-only experiment, the signal of the 
smaller, relatively slowly rotating satellite is con¬ 
densed to a smaller bandwidth that is superim¬ 
posed upon the broadband signal of the larger, of¬ 
ten rapidly rotating, primary (Fig. 1, top). Not 
all radar-observed binaries present this character¬ 
istic spectrum (e.g., where the secondary spins 
faster than the synchronous rate), but all are read¬ 
ily detected in radar images when the components 
are also resolved spatially (Fig. 1, bottom). Be¬ 
cause the spatial resolution achieved with radar 
instruments corresponds to an effective angular 
resolution of better than ~1 milliarcsecond (mas), 
there is no bias against the detection of satel¬ 
lites orbiting very close to the primary compo¬ 
nent. Multiple measurements of the range and 
frequency separations of the components over 
days of sky motion provide the geometric lever¬ 
age required to determine the orbit of the sec¬ 
ondary around the primary. This can be done 
for any orbital orientation and yields the total 
system mass, a property that is difficult to esti¬ 
mate otherwise. Other techniques involve analyz¬ 
ing spacecraft flyby and orbit trajectories (e.g.. 
Yeomans et al. 1999), measuring the Yarkovsky 
orbital drift in conjunction with thermal proper¬ 
ties (e.g., Chesley et al. 2014), or observing the 
gravitational perturbations resulting from asteroid 
encounters (e.g., Hilton 2002). 

Most binary NEA systems observed to date 
have a rapidly rotating primary and a smaller 
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Fig. 1.— Binary near-Earth asteroid (285263) 1998 
QE 2 as detected using the Arecibo planetary radar 
system. In the frequency-only spectrum showing 
echo power as a function of Doppler frequency (top), 
the narrowband echo of the tidally locked secondary 
stands out against the broadband echo of the larger, 
faster-rotating primary. In the radar image (bottom), 
the components arc spatially resolved (7.5 m/pixel). 
The vertical axis represents distance from the observer 
increasing downward. The horizontal axis is Doppler 
frequency due to the orbital and rotational motion of 
the components. Note that if one summed the pixel 
values in each column of the image, the intensity as 
a function of Doppler frequency would approximate 
the spectrum above. The secondary is roughly one- 
fourth the size of the primary (measured in the verti¬ 
cal dimension), though the Doppler breadth of the pri¬ 
mary gives the illusion of a greater size disparity. The 
shape of the secondary (inset) is distinctly nonspheri- 
cal when viewed with finer frequency resolution. 


secondary of order a few tenths the size of the 
primary (a secondary-to-primary mass ratio of 
roughly 0.001 to 0.1), whose rotation is synchro¬ 
nized to the mutual orbit period. The majority 
of primaries rotate in less than 2.8 h, though 
they range from 2.2593 h for (65803) Didy- 
mos (Pravec et al. 2006) to 4.749 h for 1998 
QE 2 (P. Pravec, pers. comm., 2013). The known 
outlier is the nearly equal-mass binary (69230) 
Hermes, whose components both appear to have 
13.894 h periods synchronized to their mutual or¬ 
bit period (Margot et al. 2006). This doubly syn¬ 
chronous configuration is most likely due to rapid 
tidal evolution (Taylor and Margot 2011). While 
the rotations of satellites in NEA binaries tend 
to be tidally locked to their orbital mean motions 
with periods typically within a factor of two of 24 
h (often resulting in the characteristic appearance 
shown in Fig. 1), about one in four radar-observed 
multiple-asteroid systems have an asynchronous 
satellite (Brozovic et al. 2011), all of which ro¬ 
tate faster than their orbital rate. Well-studied 
examples include (35107) 1991 VH (Naidu et al. 
2012), (153958) 2002 AM 31 (Taylor et al. 2013), 
(311066) 2004 DC (Taylor et al. 2008), and the 
outer satellites of both undisputed triple systems 
(153591) 2001 SN 263 (Nolan et al. 2008; Fang 
et al. 2011; Becker et al. 2015) and (136617) 
1994 CC (Brozovic et al. 2011; Fang et al. 2011). 
Of the known asynchronous satellites, all have 
wide component separations (>7 primary radii), 
translating to longer-than-typical orbital periods, 
and/or eccentric orbits (>0.05), that are either 
remnants of their formation mechanism or prod¬ 
ucts of subsequent dynamical evolution (Fang and 
Margot 2012c). 

The shortest orbital periods detected with radar 
so far are those of Didymos and 2006 GY 2 with 
P orb = 11.90lo!o2 h an d 11.7 ± 0.2 h, respec¬ 
tively (Benner et al. 2010; Brooks 2006). For 
Didymos, the semi-major axis is a = 1.181 q 02 km, 
just outside the classical fluid Roche limit of ~1 
km for equal-density components. Other sys¬ 
tems with satellites orbiting near this limit include 
2002 CE 26 and 2001 SN 263 . The significance of 
this limit is unclear, as ~100 m secondaries with 
a cohesion comparable to comet regolith or sand 
can likely survive on orbits interior to the Roche 






limit (Taylor and Margot 2010, and references 
therein). 

Inversion of a series of radar images can pro¬ 
vide a three-dimensional shape model and com¬ 
plete spin-state description given sufficient sig¬ 
nal, resolution, and orientational coverage (Hud¬ 
son 1993; Magri et al. 2007). Shape recon¬ 
struction of the larger component of (66391) 
1999 KW 4 (Ostro et al. 2006) demonstrated that 
the canonical shape of an NEA primary has 
a characteristic circular equatorial bulge, uni¬ 
formly sloped sides, and polar flattening akin to 
a spinning top. Such a shape is shared by the 
primaries of 2004 DC, 1994 CC, 2001 SN 263 , 
and (185851) 2000 DP 107 (Naidu et al. 2015), 
though some primaries have less pronounced 
equatorial belts, e.g., 2002 CE 26 and 1998 QE 2 . 
Some single asteroids have a similar shape, 
e.g., (101955) Bennu (Nolan et al. 2013) and 
(341843) 2008 EV 5 (Busch et al. 2011), but do 
not have satellites, possibly because one has not 
yet formed or has been lost in the past. Shape 
model renditions are shown in Benner et al. (this 
volume). Often the resolution of radar images 
of the smaller satellites is insufficient for shape 
inversion, but radar images suggest that the satel¬ 
lites are typically elongated, e.g., 2000 DPi 07 , 
1999 KW 4 , 2001 SN 263 , 1991 VH, and 1998 QE 2 . 

Shapes and volumes obtained from inversion 
of radar images, combined with the system mass 
derived from the orbital motion observed in radar 
images, provide the density of the system (or of 
the individual components if the mass ratio is 
measurable from reflex motion). Low densities 
of order 1 g/cm 3 (Shepard et al. 2006; Becker 
et al. 2015) to 2 g/cm 3 (Ostro et al. 2006; Bro- 
zovic et al. 2011) suggest significant internal 
macroporosity of order 50%, implying a rubble- 
pile internal structure for the components. At 
such low densities, the rapid rotation of the pri¬ 
mary places particles along the equatorial belt 
in a near-weightless environment. The combina¬ 
tion of rapid rotation, shape, and implied porosity 
and rubble-pile structure has implications for the 
formation mechanism of small multiple-asteroid 
systems (Section 4). 

While radar allows for direct, unambiguous de¬ 


tection of asteroid satellites, its range is limited. 
Because radar requires the transmission and re¬ 
ception of a signal, the strength of the received 
signal falls as the fourth power of the distance to 
the target and, thus, is best suited for detecting 
multiple-component systems passing within ~0.2 
astronomical units (au) of Earth. Satellites in the 
main asteroid belt simply tend to be too small and 
too far away to detect with present radar capabil¬ 
ities and require application of different observa¬ 
tional techniques. 

2.2 Lightcurve Observations of NEA and 
Small MBA Systems 

A photometric lightcurve is a time series of 
measurements of the total brightness of an as¬ 
teroid. Detections of binary asteroids by photo¬ 
metric lightcurve observations utilize the fact that 
the components can obscure or cast a shadow on 
one another, producing occultations or eclipses, 
respectively. The attenuations can be used to both 
reveal and characterize binaries (Fig. 2). The 
observational, analysis, and modeling techniques 
were described in Pravec et al. (2006); Scheirich 
and Pravec (2009); Scheirich et al. (2015). 

Early reports (Tedesco 1979; Cellino et al. 
1985) of asteroids suspected to be binaries on 
the basis of anomalous lightcurves (including 
(15) Eunomia, (39) Laetitia, (43) Ariadne, (44) 
Nysa, (49) Pales, (61) Danae, (63) Ausonia, (82) 
Alkmene, (171) Ophelia, and (192) Nausikaa) 
have remained largely unconfirmed despite ex¬ 
tensive follow-up searches. The first serious can¬ 
didate for detection with this technique was NEA 
(385186) 1994 AW 4 (Pravec and Hahn 1997), 
whose binary nature was confirmed by photo¬ 
metric observations in 2008 (Birlan et al. 2010). 
Since 1997, nearly 100 binaries among near-Earth 
and small main belt asteroids have been detected 
with the photometric method. The binary asteroid 
database constructed by Pravec and Harris (2007) 
(http://www.asu.cas.cz/~asteroid/binastdata.htm) 
includes data for 86 MBA and NEA binaries that 
were securely detected by photometry and for 
which basic parameters have been derived, such 
as the primary spin period, the orbital period, and 
the primary-to-secondary mean diameter ratio. A 
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Fig. 2.— Lightcurve data of (1338) Duponta, which has a secondary-to-primary diameter ratio of about 0.24. 
(a) The original data showing both lightcurve components, folded with the orbit period, (b) The orbital lightcurve 
component, derived after subtraction of the primary lightcurve component, showing the mutual events between 
components of the binary system, (c) The primary lightcurve component. Figure from Pravec et al. (2012). 


few tens of additional MBAs and NEAs are sus¬ 
pected to be binaries and await confirmation with 
more detailed observations in the future. 

Among the main findings obtained from pho¬ 
tometric observations is that binary asteroids are 
ubiquitous. They have been found among NEAs, 
Mars-crossers (MCs), and throughout the main 
belt, both among asteroids that have been identi¬ 
fied as family members and among asteroids that 
have not. Pravec et al. (2006) derived the fraction 
of binaries among NEAs larger than 300 meters 
to be 15 ± 4%. A binary fraction among MBAs 
has not been derived precisely due to less well- 
characterized observational selection effects, but 
their photometric discovery rate is similar to the 
discovery rate of binaries among NEAs. Thus, 
binaries are suspected to be as frequent among 
MBAs as they are among NEAs. There appears 


to be an upper limit on the primary diameter for 
photometrically detected binaries of about 13 km; 
the largest detected binary is (939) Isberga with 
D p = 13.4 ± 1.3 km (Carry et al. 2015). A lower 
size limit on the primary diameter D p is less clear. 
The smallest detected binary is 2000 UGn with 
D p = 0.26 ± 0.03 km (Pravec et al. 2006), but 
smaller binaries are known to exist (Section 2.1). 
Their absence in lightcurve data sets may be due 
in part to a bias against detecting small binaries in 
the initial surveys. 

Another key finding is that small binary aster¬ 
oids have, with only two or three exceptions, a 
near-critical angular momentum content (Fig. 3). 
As shown by Pravec and Harris (2007), their an¬ 
gular momentum is consistent with formation by 
fission of critically spinning parent bodies of a co¬ 
hesionless, rubble pile structure. The exceptions 
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Fig. 3.— Estimated values of the normalized total 
angular momentum content of binaries versus primary 
diameter. The quantity ol is the sum of orbital and 
spin angular momenta normalized by the angular mo¬ 
mentum of an equivalent sphere spinning at the criti¬ 
cal disruption spin rate u,i = \JiirpC /3 where p is 
the density and G is the gravitational constant. In the 
Darwin notation, = 1 corresponds to J/ J' = 0.4. 
Group A contains small NEA, MC, and MBA binaries. 
Group B consists of doubly synchronous small MBAs 
with nearly equal-size components. Group L repre¬ 
sents large MBAs with small satellites (Section 2.5). 
Two exceptional cases are the doubly synchronous as¬ 
teroids (90) Antiope and (617) Patroclus (Section 2.5). 
Figure updated from Pravec and Harris (2007). 


are the semi-wide systems (32039) 2000 JO 23 and 
(4951) Iwamoto, and possibly also (1717) Arlon, 
with orbital periods of 117 h to 360 h and super¬ 
critical total angular momentum content. 

The orbital poles of main belt binaries were 
found to have a highly anisotropic distribution, 
concentrating within 30 degrees of the poles of 
the ecliptic (Pravec et al. 2012). The preferen¬ 
tial orientations of the orbital poles suggest that 


their parent bodies or the primaries were tilted 
by the YORP effect towards the asymptotic spin 
states near obliquities 0 and 180 degrees, consis¬ 
tent with observations of single asteroids (Planus 
et al. 2011 ). 

Another significant finding is that there ap¬ 
pears to be a lower limit on the separation be¬ 
tween components of binary systems of about 
a/D p = 1.5, corresponding to an orbital period of 
11-12 h for typical densities. Lightcurve observa¬ 
tions indicate that the orbital period of Didymos 
is P orb = 11.91 ± 0.02 h (Pravec et al. 2006), 
consistent with the radar estimate. This suggests 
an orbit close to the Roche limit for strengthless 
satellites (but see prior remark about orbits inte¬ 
rior to the Roche limit). 

Photometric observations of a binary system 
over multiple apparitions can be used to detect 
a change in the separation of the components 
due to the effect on mutual event timing. An 
extensive set of photometric observations of the 
synchronous binary (175706) 1996 FG 3 obtained 
during 1996-2013 places an upper limit on the 
drift of its semi-major axis that is one order of 
magnitude less than estimated on the basis of the 
BYORP theory (Scheirich et al. 2015). This sys¬ 
tem may be in an equilibrium between BYORP 
and tidal torques as proposed for synchronous bi¬ 
nary asteroids by Jacobson and Scheeres (2011a). 

Some data sets strongly suggest the presence of 
triple asteroids. In these cases, an additional rota¬ 
tional component that does not belong to the pri¬ 
mary or the close eclipsing secondary is present 
in the lightcurve. This additional rotational com¬ 
ponent does not disappear during mutual events 
where the eclipsing close secondary is obscured 
by or in the shadow of the primary. Pravec et al. 
(2012) identified three such cases: (1830) Pog- 
son, (2006) Polonskaya, and (2577) Litva. The 
latter has been confirmed by direct imaging ob¬ 
servations of the third body (second satellite) on 
a wide orbit (Merline et al. 2013). 

Other data sets reveal the existence of paired 
binaries/triples. Two such cases have been pub¬ 
lished: the pair composed of (3749) Balam and 
2009 BR 60 (Vokrouhlicky 2009, and references 
therein) and the pair composed of (8306) Shoko 
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and 2011 SR 158 (Pravec et al. 2013). Balam 
is a confirmed triple, with a distant satellite de¬ 
tected by direct imaging (Merline et al. 2002a) 
and a close satellite detected by lightcurve obser¬ 
vations (Marchis et al. 2008d). Shoko is a sus¬ 
pected triple as well: Using lightcurve observa¬ 
tions, Pravec et al. (2013) detected an eclipsing, 
synchronous close satellite with P or b = 36.2 h 
and a third rotational component attributed to an 
outer satellite. 

While the population of binary NEAs and small 
MBAs is composed primarily of synchronous 
systems, and secondarily of asynchronous sys¬ 
tems with low secondary-to-primary size ratios 
(. D s /D p < 0.5), doubly synchronous binaries 
with nearly equal-size components also exist 
(Fig. 4). Nine such systems with D s jD v > 0.7 
and orbital periods between 15 h and 118 h have 
been reliably identified in the main belt (e.g., 
Behrend et al. 2006; Kryszczynska et al. 2009, 
see also the Pravec and Harris binary database 
described above). 

Another important observation is that, with the 
exception of doubly synchronous systems, all bi¬ 
naries have unelongated, near-spheroidal primary 
shapes, as evidenced by their low primary am¬ 
plitudes not exceeding 0.3 mag (when corrected 
to zero phase angle). This suggests that their 
primaries may have shapes similar to the top¬ 
like shapes that have been observed for 1999 
KW 4 (Ostro et al. 2006) and several other bina¬ 
ries by radar. 

All the properties revealed by photometric ob¬ 
servations indicate that binary systems among 
NEAs and small MBAs were formed from crit¬ 
ically spinning cohesionless parent bodies, with 
YORP as the predominant spin-up mechanism. 
This finding is consistent with the fact that the 
observed 0.2-13 km size range of binaries cor¬ 
responds to the size range where the spin barrier 
against asteroid rotations faster than about 2.2 h 
has been observed (e.g., Pravec et al. 2007). 

Although lightcurve observations provide pow¬ 
erful constraints on binaries, there are limita¬ 
tions. Detection of mutual events requires an 
edge-on geometry and observations at the time of 
the events, such that some binaries remain unde- 
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Fig. 4.— Primary rotation period versus primary di¬ 
ameter. Groups A, B, and L arc defined in the caption 
of Fig. 3. Three doubly synchronous asteroids with 
nearly equal-size components lie isolated in the plot: 
(69230) Hermes on the left and (90) Antiope and (617) 
Patroclus on the right of group B. Note that members 
of group A cluster near the disruption spin limit for 
strengthless bodies. Figure from Pravec and Harris 
(2007). 

tected (e.g., (69230) Hermes during its 2003 ap¬ 
parition). Small satellites also escape detection 
because their effect on the lightcurve is not mea¬ 
surable (e.g., satellites with D s /D p < 0.17 re¬ 
main undetected if the minimum detectable rel¬ 
ative brightness attenuation is ~0.03 mag). The 
probability of mutual event detection is larger at 
smaller semi-major axes (expressed in units of 
primary radius) and at larger size ratios, result¬ 
ing in observational biases (e.g., Pravec et al. 
2012). Finally, lightcurve observations yield rel¬ 
ative, not absolute, measurements of orbital sepa¬ 
rations. Detection of small or distant secondaries 
and direct measurement of orbital separation must 
instead rely on other observational techniques. 

2.3 Lightcurve Observations of Asteroid Pairs 

Vokrouhlicky and Nesvomy (2008) reported 
evidence for pairs of MBAs with bodies in each 
pair having nearly identical heliocentric orbits. 
Because chance associations can be ruled out, 
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the asteroids in each pair must be genetically re¬ 
lated. Quantifying the difference in orbital pa¬ 
rameters is accomplished with a metric d that 
corresponds roughly to the relative velocity be¬ 
tween the bodies at close encounter. Vokrouh¬ 
licky and Nesvorny (2008) identified 44 asteroid 
pairs (excluding family members) with a distance 
between the orbits of their components amount¬ 
ing to d < 10 m/s. They showed that, when in¬ 
tegrated backwards in time, the orbits converge at 
a certain moment in the past with a physical dis¬ 
tance much less than the radius of the Hill sphere 
and with a low relative velocity on the order of 
1 m/s. 

Pravec and Vokrouhlicky (2009) developed a 
method to identify probable asteroid pairs by se¬ 
lecting candidate pairs with a similar distance cri¬ 
terion, then computing the probability that each 
candidate pair emerged as a result of a coinci¬ 
dence between two unrelated asteroids. They 
identified 72 probable asteroid pairs, reproducing 
most of the 44 previously known pairs. Most of 
the new candidates were later confirmed to be real 
pairs using backward integrations of their helio¬ 
centric orbits. 

Vokrouhlicky and Nesvorny (2008) proposed 
a few possible formation mechanisms for the as¬ 
teroid pairs: collisional disruption, rotational fis¬ 
sion, and splitting of unstable asteroid binaries. 
Pravec et al. (2010) conducted a survey of the 
rotational properties of asteroid pairs, and they 
found a strong correlation between the primary 
rotational periods and the secondary-to-primary 
mass ratio (Fig. 5). They showed that this correla¬ 
tion fits precisely with the predictions of a model 
by Scheeres (2007) in which a parent body with 
zero tensile strength undergoes rotational fission. 
The model predicts that primaries of low mass ra¬ 
tio pairs (q < 0.05) have not had their spin sub¬ 
stantially slowed down in the separation process 
and should rotate rapidly with frequencies close 
to the fission spin rate. The observed periods are 
between 2.4 and 5 h. Primaries of medium mass 
ratio pairs (q = 0.05 to ~ 0.2) have had their 
spin slowed down according to the model because 
a substantial amount of angular momentum was 
taken away by the escaped secondary. This trend 
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Fig. 5.— Primary rotation periods versus mass ratios 
of asteroid pairs. The mass ratio values were estimated 
from the differences between the absolute magnitudes 
of the pair components, AH. Circles arc data points 
with quality code rating U p = 3, meaning a precise 
period determination. Diamonds arc data points with 
Up = 2, which arc somewhat less certain estimates. 
Error bars arc one standard deviation. The data match 
the predictions (curves) of a model of rotational fission 
with a few adjustable parameters. In the model, Tj n ; is 
the binary system's initial orbit semi-major axis, q/_ is 
the normalized total angular momentum of the system 
(Fig. 3), and a p , b p , c p arc the long, intermediate, and 
short axis of the dynamically equivalent equal mass 
ellipsoid of the primary. All models shown assume 
b p /cp = 1.2. The dashed curve shows the best-fit 
model with ol = 1.0, a p /b p = 1.4 and A- m \/b p = 3. 
Solid curves represent upper and lower limiting cases 
with ol = 0.7—1.2. Figure updated from Pravec et al. 
( 2010 ). 

is observed in the data (Fig. 5). Finally, high mass 
ratio pairs with q > 0.2 should not exist, as the 
free energy in the proto-binary system formed by 
rotational fission would be negative and the com¬ 
ponents would be unable to separate. Observa¬ 
tions mostly corroborate this prediction: all 32 
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pairs in the sample of Pravec et al. (2010) were 
found to have a mass ratio < 0.2. However, an ex¬ 
panded photometric survey with 64 asteroid pairs 
observed between 2012 and the date of this writ¬ 
ing reveals 3 pairs with high mass ratio (q > 0.5). 
Their formation requires an additional supply of 
angular momentum. Another important finding 
by Pravec et al. (2010) is that the primaries of as¬ 
teroid pairs have lightcurve amplitudes that im¬ 
ply shapes with a broad range of elongations, i.e., 
unlike the primaries of binaries (Sections 2.1 and 
2.2), the primaries of asteroid pairs do not tend to 
be nearly spheroidal. 

2.4 Spectral Observations of Asteroid Pairs 

Colorimetric and spectral observations of about 
20 asteroid pairs indicate that members of an as¬ 
teroid pair generally have similar spectra (Duddy 
et al. 2012; Moskovitz 2012; Duddy et al. 2013; 
Polishook et al. 2014a; Wolters et al. 2014). In 
some pairs, the authors observed subtle spec¬ 
tral differences between the components and at¬ 
tributed them to a larger amount of weathered ma¬ 
terial on the surface of the primary. In two pairs, 
they observed somewhat more significant spec¬ 
tral differences. For the pair (17198)-(229056), 
both Duddy et al. (2013) and Wolters et al. (2014) 
found that the primary is redder, i.e., it has a 
somewhat higher spectral slope than the sec¬ 
ondary in the observed spectral range 0.5-0.9 /im. 
It is unclear why their spectra differ despite a 
strong dynamical link between the two aster¬ 
oids. For the pair (19289)-(278067), Wolters 
et al. (2014) observed a spectral difference simi¬ 
lar to that seen in (17198)—(229056), but Duddy 
et al. (2013) observed very similar spectra. Cross- 
validation of the methods or additional observa¬ 
tions, perhaps rotationally resolved, are needed to 
resolve the discrepancy. 

2.5 Direct Imaging of MBA and Trojan Sys¬ 
tems 

Direct imaging of asteroids can reveal the pres¬ 
ence of satellites and, following the long tradition 
of orbit determination of binary stars and plan¬ 
etary satellites, lead to estimates of orbital pa¬ 
rameters (Fig. 6). This observing mode remains 


challenging because the satellites are generally 
much smaller and fainter than their respective pri¬ 
maries and because most satellites known to date 
orbit at angular separations below 1 arcsecond. 
Satellite discoveries have therefore followed the 
development of adaptive optics (AO), and recent 
advances have enabled the detection of asteroid 
satellites that had remained undetected in prior 
searches. 

Instruments must have sufficient contrast and 
resolving power to detect asteroid satellites with 
direct imaging. For a 50-100 km diameter aster¬ 
oid in the main belt orbited by a satellite a few 
km across, the typical angular separation is gen¬ 
erally less than an arcsecond with a contrast of 5 
to 10 magnitudes (computed as 2.5 log (F p /F s ), 
where F is the flux and p and s indicate primary 
and secondary, respectively). 

In some situations, direct images can actually 
resolve the primary. A 50-100 km diameter aster¬ 
oid at 2 au subtends 34-68 mas while the diffrac¬ 
tion limit of a 10 m telescope at a typical imag¬ 
ing wavelength of 1.2 pm is about 30 mas. Al¬ 
though the diffraction limit is not reached, it can 
be approached with high-performance AO instru¬ 
ments in excellent conditions. With a sequence of 
disk-resolved images that provide sufficient ori¬ 
entational coverage, it is possible to estimate the 
3D shape of the primary. This enables volume 
and density determinations. 

Instruments capable of meeting the contrast 
and resolution requirements include the Hub¬ 
ble Space Telescope (HST) and large (10 m 
class) ground-based telescopes equipped with 
AO. Spacecraft encounters provide an opportu¬ 
nity to detect small satellites at small separations 
because of proximity to the target and the ab¬ 
sence of the point spread function halo that affects 
ground-based AO instruments. 

At the time Asteroids III was published, MBA 
satellite discoveries included one by spacecraft 
((243) Ida), one by HST ((107) Camilla), and 6 
by ground-based AO instruments. Since then, 
ground-based AO instruments have been respon¬ 
sible for almost all large MBA satellite discov¬ 
eries: (121) Hermione (Merline et al. 2002b), 
(379) Huenna (Margot 2003), (130) Elektra 
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Fig. 6.— Satellite detection by direct imaging with adaptive optics (AO), (a) Image of asteroid (41) Daphne 
(Vmag=10) obtained with a ground-based AO-fed camera (NACO at ESO VLT, 5 s exposure), (b) Same image 
after subtraction of the flux from the primary, enabling more accurate measurements of the flux and position of 
the secondary, (c) Orbit determination. The relative positions of the satellite from VLT/NACO and Keck/NIRC2 
images arc indicated. Figure adapted from Carry (2009). 


(Merline et al. 2003c), a second satellite to (87) 
Sylvia (Marchis et al. 2005b) and to (45) Euge¬ 
nia (Marchis et al. 2007), (702) Alauda (Rojo and 
Margot 2007), (41) Daphne (Conrad et al. 2008), 
two satellites to (216) Kleopatra (Marchis et al. 
2008b) and (93) Minerva (Marchis et al. 2009), 
and (317) Roxane (Merline et al. 2009). The wide 
binaries (1509) Esclangona (Merline et al. 2003a) 
and (4674) Pauling (Merline et al. 2004), which 
are small asteroids in our classification, have also 
been identified using AO-fed cameras. HST en¬ 
abled detections of two additional wide binaries: 
(22899) 1999 TO M (Merline et al. 2003b) and 
(17246) 2000 GL 74 (Tamblyn et al. 2004), both of 
which are small MBAs. No satellites have been 
discovered around any of the 7 asteroids recently 
visited by spacecraft: (4) Vesta, (21) Lutetia, 
(2867) Steins, (4179) Toutatis, (5535) Annefrank, 
(25143) Itokawa, and (132524) APL. The num¬ 
ber of known large MBAs with satellites is now 
16, which includes the only known large dou¬ 
bly synchronous system, (90) Antiope (Merline 
et al. 2000; Michalowski et al. 2004; Descamps 
et al. 2007, 2009). The fraction of large MBAs 
with satellites is difficult to estimate because of 
a complex dependence of satellite detectability 
on primary-to-secondary angular separation and 
primary-to-secondary flux ratio. However, be¬ 
cause several independent programs have sur¬ 


veyed over 300 large MBAs, it is likely that the 
abundance of binaries in large MBAs is substan¬ 
tially smaller than the ~16% abundance in NEAs 
and small MBAs. 

Properties of large MBA binaries and triples 
are summarized in Figs. 7 and 8. With the excep¬ 
tion of the nearly equal-mass binary (90) Antiope, 
the known satellites have secondary-to-primary 
mass ratios between 10~ 6 and 10~ 2 . All have 
orbital periods between 1 and 5.5 days, except 
(379) Huenna, whose orbit has a period of ~88 
days and an eccentricity of ~0.2 (Marchis et al. 
2008c). Many orbits have near-zero eccentric¬ 
ity (e.g., Marchis et al. 2008a), likely the result 
of tidal damping, but the inner satellites of triples 
generally have non-zero eccentricities. These 
eccentricities may have originated when orbits 
crossed mean motion resonances while tidally ex¬ 
panding (e.g., Fang et al. 2012). 

At first glance, large MBA densities appear to 
cluster in two groups, between 1 and 2 g/cm 3 and 
above 3 g/cm 3 . However, interpretations are lim¬ 
ited by the possibility of systematic errors, in¬ 
cluding overestimates of volumes and underesti¬ 
mates of densities (Pravec and Harris 2007). Be¬ 
cause volume uncertainties almost always dom¬ 
inate the error budget for binary asteroid densi¬ 
ties (e.g., Merline et al. 2002c; Carry 2012), it is 
important to assess the realism of uncertainties 
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Fig. 7.— Properties of large MBA binaries and 

triples, excluding the doubly synchronous (90) An- Fig. 8.—Properties of satellites of large MBAs, ex- 
tiope. Error bars or upper limits, when available, are eluding outliers (90) Antiope and (379) Huenna (see 
shown. Figures 7 and 8 are based on data compiled text ). Satellites of (243) Ida and (317) Roxane, whose 
by Johnston (2014) from references therein. orbits arc not well known, arc not shown. 
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associated with volume determinations. Some 
published density values should be regarded with 
caution because overconfidence in the fractional 
uncertainty of volume estimates has led to un¬ 
derestimates of bulk density uncertainties. The 
platinum standard of an orbiting spacecraft yields 
densities with ~1% accuracy. The gold standard 
of radar observations where tens of images with 
hundreds or thousands of pixels per image are 
used to reconstruct a detailed 3D shape model 
yields volumes (and densities) with ~10% accu¬ 
racy. In contrast, AO images contain at most a few 
independent resolution cells of the target aster¬ 
oid. Shape reconstructions based on AO images 
and/or lightcurve data may not routinely yield 
volume accuracies at the 10% level, although one 
analysis reached that level (Carry et al. 2012). In 
the absence of precise volume information, one 
might be tempted to infer bulk densities from 
the theory of fluid equilibrium shapes, but this 
approach is problematic (Holsapple 2007; Harris 
et al. 2009). 

In the Jupiter trojan population, one satellite 
to (624) Hektor has been reported (Marchis et al. 
2006b) since the discovery of the first trojan satel¬ 
lite to (617) Patroclus (Merline et al. 2001). These 
are the only trojans confirmed to have satellites 
in spite of several active search programs. The 
apparent low abundance of binary trojans is in¬ 
triguing and, if confirmed, may provide addi¬ 
tional support for the idea that Jupiter trojans 
originated in the trans-Neptunian region (Mor- 
bidelli et al. 2005; Levison et al. 2009) where 
they experienced a different collisional environ¬ 
ment than in the main belt of asteroids. (624) 
Hektor has a satellite in a ~3-day orbit that is ec¬ 
centric (~0.3) and inclined (~50°) with respect to 
Hektor’s equator (Marchis et al. 2014). (617) Pa¬ 
troclus is unusual because it has two components 
of similar size in a relatively tight (~680 km) or¬ 
bit, with a normalized total angular momentum 
exceeding that available from fission of a single 
parent body (Marchis et al. 2006a). 

In the trans-Neptunian region, 14 and 64 binary 
systems have been discovered with AO and HST, 
respectively (Johnston 2014). The apparent larger 
abundance of binary TNOs in the cold classical 


belt may be due to a different dynamical environ¬ 
ment and formation mechanism (Section 5). 

Objects in the trojan and TNO populations are 
generally too faint for AO observations in natural 
guide star (NGS) mode, in which the science tar¬ 
get is also used to measure the properties of the 
wavefront and command the deformable mirror. 
These objects can be observed in appulse when 
their sky position happens to be within < 1 ar- 
cminute of a bright star. The advent of laser guide 
star (LGS) adaptive optics has been an impor¬ 
tant development that has freed the observer from 
finding such chance alignments and has opened 
up a larger fraction of the sky for observation of 
faint objects. Even with LGS, however, the avail¬ 
ability of a tip-tilt star (Rmag < 18) within < 1 
arcminute of the target is still required. 

High-resolution and high-contrast imaging 
capabilities are aggressively sought by instru¬ 
ment builders, in part to enable direct imaging 
of exoplanets. Cameras equipped with high- 
performance AO are currently being installed 
or commissioned on large ground-based tele¬ 
scopes: HiCIAO on Subaru, GPI on Gemini, and 
SPHERE at the ESO VLT. These instruments will 
improve the ability to detect faint satellites orbit¬ 
ing close to their respective primaries. However, 
in most cases, asteroids fall in the faint-end range 
of these instrument capabilities. The next gener¬ 
ation of large telescopes (~30 m diameter) such 
as the Thirty Meter Telescope (TMT) and Eu¬ 
ropean Extremely Large Telescope (E-ELT) wifi 
provide an improvement in sensitivity by a factor 
of ~10 and in angular resolution by a factor of 
~3 compared to current 10 m telescopes. With 
the anticipated development of AO capabilities at 
shorter wavelengths, the second generation of in¬ 
struments at these facilities is expected to provide 
improvements in angular resolution by a factor of 
~5. Such instruments may allow detection of the 
small MBA binaries that are currently beyond the 
reach of direct imaging instruments. In many of 
these systems, the components are separated by 
only a few mas and the size ratios are larger than 
in large MBA binaries, resulting in flux ratios 
closer to unity. 
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2.6 Spectral Observations of MBA and Trojan 
Systems 

It is generally difficult to separate the light 
emitted or reflected from the secondary from that 
of the brighter primary. Nevertheless, such ob¬ 
servations can be attempted when the secondary 
happens to be at a large angular separation from 
the primary, when the system is undergoing mu¬ 
tual events, or with the help of an integral field 
spectrograph. 

Spectra of (22) Kalliope and its satellite Li¬ 
nus in the 1-2.4 /im region appear to be sim¬ 
ilar (Laver et al. 2009), which the authors at¬ 
tribute to satellite formation after a major impact 
on the precursor body. Observations of both com¬ 
ponents of (90) Antiope in the same spectral re¬ 
gion also shows surface reflectances that are sim¬ 
ilar (Marchis et al. 2011). The spectrum of (379) 
Huenna is characteristic of C-type asteroids and 
the secondary does not exhibit a significantly dif¬ 
ferent taxonomic type (DeMeo et al. 2011). Both 
components of (809) Lundia are consistent with a 
V-type classification (Birlan et al. 2014). 

In the mid-infrared, Spitzer observations of the 
trojan (617) Patroclus, including during mutual 
events, provided size estimates for its components 
and a thermal inertia of 20 ± 15 J s -1 / 2 Kr 1 m~ 2 
(Mueller et al. 2010). Spitzer observations com¬ 
bined with photometric results in the visible 
yielded size and albedo estimates for (624) Hek- 
tor (Emery et al. 2006). Spitzer observations of 
these and other binaries did not resolve the bina¬ 
ries and results typically cannot be compared to 
observations that place many resolution elements 
on individual components. One exception is 
2000 DPi 07 , where analysis of Spitzer data yields 
a system density of 0.9 ± 0.3 g/cm 3 (Marchis 
et al. 2012) and the radar results indicate 1.4 ± 
0.2 g/cm 3 (Naidu et al. 2015). 

2.7 Stellar Occultations of MBA and Trojan 
Systems 

Stellar occultations provide a way of detect¬ 
ing components of a multiple-asteroid system, of 
placing bounds on component sizes, and of ob¬ 
taining the relative positions of components on 


the plane of the sky. A recording of star light as 
a function of time shows a deep extinction when 
a target body crosses the line of sight between the 
observer and the star. This can be interpreted in 
terms of a chord on the apparent disk of the tar¬ 
get body projected on the plane of the sky. If 
several observers are placed across the occupa¬ 
tion path on the surface of the Earth, multiple 
chords can be obtained, and the size and shape 
of the target projected on the sky can be recon¬ 
structed (Fig. 9). When two or more components 
are present, it is also possible to measure their rel¬ 
ative position. While the reliability of this tech¬ 
nique was disputed a decade ago due to the lack 
of digital recordings, the availability of low-cost 
cameras and global positioning systems has en¬ 
abled a dramatic improvement in the precision of 
timing reports. Stellar occultations have become 
an important observational tool for the study of 
binary asteroids. 

Early reports (e.g., Binzel and van Flandem 
1979) of asteroids suspected to be binaries on the 
basis of occultation data (including (3) Juno, (6) 
Hebe, (9) Metis, (12) Victoria, (18) Melpomene, 
(146) Lucina, and (532) Herculina) have re¬ 
mained largely unconfirmed despite extensive 
follow-up searches. However, it is likely that 
the outer satellite of (216) Kleopatra was de¬ 
tected during a 1980 occultation (Dunham 1981; 
Descamps et al. 2011). The detection of a satel¬ 
lite around the trojan (911) Agamemnon has been 
suggested (Timerson et al. 2013) but not yet con¬ 
firmed. The occultation technique has also been 
used to detect rings around the centaur (10199) 
Chariklo (Braga-Ribas et al. 2014). 

One strength of the stellar occultation tech¬ 
nique lies in the fact that the observability of the 
event depends mainly on the brightness of the star 
and not of the asteroid or satellite. Stellar occul¬ 
tations can thus be used to detect small (km size) 
satellites, even those that are close to the primary 
and that would remain undetected in direct imag¬ 
ing. 

Another strength of the technique is the po¬ 
tential for high-precision measurements. Stellar 
occultations are based on time-series photometry. 
Given a sufficiently high cadence (e.g., 10-30 im- 
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ages per second), it is possible to obtain a preci¬ 
sion of a few mas on the relative position of bi¬ 
nary components, which is 5 to 10 times better 
than with direct imaging with current instrumen¬ 
tation. 

Finally, well-sampled stellar occultations al¬ 
low for recovery of the size and apparent shape 
of asteroids and their satellites, whereas optical 
lightcurves and direct imaging observations pro¬ 
vide primarily the diameter ratio of the compo¬ 
nents and more limited shape information. So far, 
four successful observations of satellite size and 
shape have been reported: Linus, satellite of (22) 
Kalliope (Descamps et al. 2008), Romulus, the 
outer satellite of (87) Sylvia (Berthier et al. 2014), 
and both components of the equal-size binaries 
(90) Antiope (Bartczak et al. 2014) and (617) Pa- 
troclus (Buie et al. 2014). 



X (km) 


Fig. 9.— The apparent shape of Linus, a satellite 
of (22) Kalliope (Margot and Brown 2003), detected 
by stellar occultations. In this analysis, the profile of 
the satellite (solid curve) fitted to the observed chords 
(straight lines) yields an equivalent diameter of 30 ± 
6 km. Dashed curves show the corresponding uncer¬ 
tainty of the fitted profile, and dashed lines show neg¬ 
ative detections. Figure adapted from Descamps et al. 
(2008). 

Despite all of these strengths, there remains a 
relatively low number of well-covered stellar oc- 
cultation events. This is due, in part, to the re¬ 
quirement of successful observations at many sta¬ 
tions. Owing to uncertainties on both the star 
and asteroid positions, the occultation path can 


shift by several tens or even hundreds of km on 
Earth compared to the prediction. Observers must 
therefore spread geographically to cover an event, 
but the detection of a satellite by several stations 
requires a fine grid of observers. 

The situation is, however, expected to improve 
dramatically with the availability of the Gaia stel¬ 
lar catalog and better asteroid orbits (Tanga and 
Delbo 2007). Predictions of the occultation paths 
(for the center of mass) will be accurate to a few 
km, and the main source of uncertainty will be¬ 
come the prediction of the relative position of the 
satellite around the primary. 

2.8 Other Observations 

There have been several attempts to use ground- 
based interferometers to measure the angular sep¬ 
aration of binary systems (Delbo et al. 2009; 
Carry et al. 2015). However, asteroid satellites 
are too faint for current interferometers operating 
in the visible and near-infrared and at the edge of 
detection in the mid-infrared. Future instrumen¬ 
tation may allow such observations. There are 
also prospects for observations with the ALMA 
sub-millimeter array (Busch 2009). 

3. DYNAMICS 

In parallel with advances in instrumentation 
and observing capabilities, the field has seen 
tremendous developments in understanding the 
dynamical processes that affect asteroid systems. 
This has been enabled in large part by the avail¬ 
ability of detailed shape models and orbital pa¬ 
rameters, by the need to model the dynamics of 
newly discovered triple systems, and by the desire 
to understand formation and evolution processes. 

A non-exhaustive list of some dynamical prob¬ 
lems that have been explored since Asteroids III 
includes the stability of asteroid satellite orbits 
(Scheeres 2002; Frouard and Compere 2012), 
the dynamics around triaxial bodies (Scheeres 
2009a), the fate of asteroid ejecta (Scheeres 
2007), the formation of contact binaries via dy¬ 
namical evolution (Scheeres 2009a; Taylor and 
Margot 2011, 2014), the genesis of eccentric and 
mutually inclined orbits (Fang et al. 2011; Fang 
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and Margot 2012c), the orbital determination of 
triple systems using point-mass approximations 
(Marchis et al. 2010) and full N-body calcula¬ 
tions (Fang et al. 2012), the influence of Kozai 
cycles on binaries (Perets and Naoz 2009; Fang 
and Margot 2012b), the effects of close planetary 
encounters on mutual orbits (Fang and Margot 
2012a) and spin states (Takahashi et al. 2013), 
the complex spin-orbit interactions with irreg¬ 
ular component shapes (Scheeres et al. 2006), 
including the libration and irregular rotation of 
secondaries (Naidu and Margot 2015), the in¬ 
fluence of internal structure (Goldreich and Sari 
2009), material properties (Taylor and Margot 
2011) and nonspherical shapes (Taylor and Mar¬ 
got 2014) on tidal evolution, the possibility of 
tidal saltation (Harris et al. 2009; Fahnestock 
and Scheeres 2009), the possibility of significant 
radiative evolution (Cuk and Burns 2005; Cuk 
2007; Cuk and Nesvorny 2010; McMahon and 
Scheeres 2010a,b), and the possibility of a stable 
equilibrium between tidal and radiative evolution 
(Jacobson and Scheeres 2011a). 

Several radar data sets provide exquisite con¬ 
straints for dynamical studies. Reflex motion 
has been measured for 2000 DP 107 (Margot et al. 
2002; Naidu et al. 2015), 1999 KW 4 (Ostro et al. 
2006), and 1991 VH (Naidu et al. 2012), allow¬ 
ing masses of individual components to be de¬ 
termined. Because detailed component shapes 
are also available, one can fully model the sys¬ 
tem dynamics and study spin-orbit coupling in 
detail (Scheeres et al. 2006; Fahnestock and 
Scheeres 2008; Naidu and Margot 2015). One 
finding from this work is that even moderately 
elongated secondaries on mildly eccentric or¬ 
bits are likely to experience chaotic rotation that 
substantially affect binary evolution timescales 
(Fig. 10). 

4. SMALL ASTEROIDS: SYNTHESIS 
4.1. Rotational Fission Hypothesis 

With the exception of the doubly synchronous 
binary asteroid systems, the primary asteroids 
of all small binary systems are rapidly rotating 



Fig. 10.— Surface of section plot showing the pos¬ 
sible rotational regimes of the ~200 m secondary of 
1991 VH (secondary elongation a/b = 1.5 and mu¬ 
tual orbit eccentricity e = 0.05). The plot shows the 
angle between the long axis and the line of apsides 
of the mutual orbit, 9 P , against its time derivative, 
Op, normalized by the mean motion, n, at each peri- 
center passage. Five trajectories arc illustrated (from 
top to bottom: non-resonant quasi-periodic, periodic, 
chaotic, periodic, periodic). While trapped in the sea 
of chaos, the secondary experiences torques on its per¬ 
manent deformation that result in a highly variable 
spin rate, preventing BYORP-type evolution. Figure 
from Naidu and Margot (2015). 

(within a factor of only a few of the critical dis¬ 
ruption spin limit for bodies with no shear or ten¬ 
sile strength u>d = \f An pG/?>). Furthermore, al¬ 
most all known small binary asteroids have high 
angular momentum contents (Pravec and Har¬ 
ris 2007). These characteristics are not consis¬ 
tent with formation following a sub-catastrophic 
impact, capture through a three-body interaction 
in the near-Earth or main belt, or capture af¬ 
ter a catastrophic impact. Instead, they are in¬ 
dicative of formation from a rotational fission 
event (e.g., Margot et al. 2002; Pravec and Har¬ 
ris 2007). The rotational fission hypothesis posits 
that a parent asteroid can be torqued to a rota¬ 
tion rate so great that the centrifugal accelerations 
overpower the gravitational accelerations holding 
a strengthless asteroid together (Weidenschilling 
1980). It is possible that some small asteroids 
have cohesive or molecular strength in addition 
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to self-gravity (e.g., Rozitis et al. 2014). In these 
cases, the centrifugal accelerations must over¬ 
come these additional forces in order for the as¬ 
teroid to fission (Pravec and Harris 2000; Sanchez 
and Scheeres 2014). At rapid rotation rates, loose 
surface material can flow from high-latitude re¬ 
gions to the equator along potential gradients (Os- 
tro et al. 2006). It has been shown that rota¬ 
tional acceleration could trigger local slope fail¬ 
ures and landslides, which can form the canon¬ 
ical top shape and equatorial bulge seen on pri¬ 
mary components in small multiple-asteroid sys¬ 
tems (Walsh et al. 2008; Harris et al. 2009). 

Bottke et al. (2002) proposed a YORP-induced 
rotational fission hypothesis. It has since been 
shown that the YORP effect controls the ro¬ 
tational acceleration of small asteroids (Bottke 
et al. 2006; Marzari et al. 2011) and naturally 
explains the period distribution among small as¬ 
teroids (Pravec et al. 2008; Rossi et al. 2009; Pol- 
ishookand Brosch 2009). Furthermore, including 
the YORP-induced rotational fission hypothesis 
in size-frequency distribution models improves 
the agreement with observations (Jacobson et al. 
2014a). The observed characteristics of the sys¬ 
tems described in Sections 2.1-2.3 as well as ther¬ 
mal inertia observations (Delbo et al. 2011) are 
consistent with a binary formation mechanism 
that involves spin-up and mass shedding. The 
YORP-induced rotational fission hypothesis is the 
leading candidate for explaining the formation of 
binaries, triples, and pairs among small asteroids. 

4.2. Asteroid Pairs 

The YORP effect can increase the spin rate of 
asteroids beyond the critical disruption spin limit, 
thereby triggering rotational fission. In actuality, 
there is some uncertainty regarding the spin rate 
at which disruption occurs—there may be failure 
and deformation before fission (Walsh et al. 2008; 
Sanchez and Scheeres 2011; Cotto-Figueroa et al. 
2013). The critical disruption spin limit also de¬ 
pends on the detailed shapes, masses, interlock¬ 
ing nature of the interior components and any 
cohesive forces (Scheeres 2007, 2009b; Sanchez 
and Scheeres 2014). Despite ignoring these de¬ 
tails, simple calculations provide a rotational fis¬ 


sion model that can be compared directly and suc¬ 
cessfully with observations. 

If a spherical approximation of each compo¬ 
nent is made, then the rotational breakup spin 
rate u q necessary for fission as a function of the 
secondary-to-primary mass ratio q is (Scheeres 
2007): 

= (1> 

This is the exact solution for two spheres resting 
on each other with a mass ratio of q and rotating 
about the axis of maximum moment of inertia. 

The spherical component model described 
above demonstrates the important reality that the 
larger the mass ratio q of the two future binary 
members, the slower the required rotation rate 
necessary to create the binary system. This slower 
required rotation rate translates into a small initial 
free energy for the ensuing binary system. The 
free energy Ef is the energy that is accessible to 
the different energy reservoirs in the system, in¬ 
cluding the rotation states of each member and 
the orbit. It does not include the internal bind¬ 
ing energy of each object. The free energy is 
an important quantity because it determines the 
boundedness of the system. Bound systems have 
negative free energy, while unbound systems have 
positive free energy. An unbound binary system 
implies that the system is capable of disruption 
but does not imply that the system will disrupt. 
For the idealized case of two spheres, the free 
energy can be expressed as (Scheeres 2007): 


E f - 


2npu 2 d Rl 

15 


/(<?), 


( 2 ) 


where R p is the radius of the primary and f(q) is 
an algebraic, monotonically decreasing function 
for 0 < q < 1. For the equation above corre¬ 
sponding to two spheres, the function crosses zero 
when q & 0.204. Similar equations can be written 
for any two component shapes, but q ~ 0.2 re¬ 
mains near the binding energy transition point, so 
the model uses this point as a simple approxima¬ 
tion. This crossing point divides bound systems 
with negative energy and mass ratios q > 0.2 
and unbound systems with positive energy and 
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mass ratios q < 0.2. Because of this fundamen¬ 
tal difference, high mass ratio q > 0.2 and low 
mass ratio q < 0.2 binary systems evolve dif¬ 
ferently (Scheeres 2009a; Jacobson and Scheeres 
2011b). Primarily, positive energy low mass ra¬ 
tio systems will chaotically explore orbital phase 
space until the majority find a disruption trajec¬ 
tory creating an asteroid pair; this evolutionary 
route is unavailable to high mass ratio systems. 

The asteroid pair population provides a natu¬ 
ral laboratory to test this relationship (Scheeres 
2007; Vokrouhlicky and Nesvorny 2008). Pravec 
et al. (2010) examined many asteroid pair sys¬ 
tems and measured the rotation rate of the pri¬ 
mary and the absolute magnitude difference be¬ 
tween the pair members. These two quantities 
should follow a simple relationship related to u q , 
although many of the ignored details mentioned 
at the beginning of this section can move aster¬ 
oids away from this relationship. Indeed, Pravec 
et al. (2010) discovered that asteroid pairs do fol¬ 
low this relationship (Fig. 5). Furthermore, they 
found that the large members of asteroid pairs 
have a broader range of elongations than the pri¬ 
maries of binary systems, consistent with the find¬ 
ings of Jacobson and Scheeres (2011b) that pro¬ 
late primaries are less likely to remain in a bound 
binary system after rotational fission. Thus, there 
is strong evidence to support the hypothesis that 
asteroid pairs are the products of rotational fis¬ 
sion. 

Asteroid pairs continue to be a fertile obser¬ 
vational landscape. Since dynamical integrations 
can derive the “birthdate” of such systems, ob¬ 
servers can test ideas regarding space weathering 
timescales and YORP evolution after fission (Pol- 
ishook et al. 2014a; Polishook 2014). Along with 
binary systems, the surfaces of asteroid pairs may 
provide clues in the future regarding the violence 
of the rotational fission process (Polishook et al. 
2014b). 

4.3. Binary and Triple Systems 

Jacobson and Scheeres (201 lb) showed that af¬ 
ter rotational fission there are a number of possi¬ 
ble outcomes. Their numerical studies produced 
the evolutionary flow chart shown in Fig. 11; 


many of these outcomes were also found by Fang 
and Margot (2012c). The high and low mass ra¬ 
tio distinction for rotational fission emphasized 
above plays an important role in distinguishing 
the two evolutionary pathways. Along the high 
mass ratio pathway, both binary members tidally 
synchronize and then evolve according to the BY- 
ORP effect. 

Along the low mass ratio pathway, the bi¬ 
nary system is unbound. Since these systems are 
chaotic, many are disrupted and become aster¬ 
oid pairs. During this chaotic binary state, the 
secondary can often go through rotational fission 
itself, although this rotational fission is torqued 
by spin-orbit coupling (Fig. 10) rather than the 
YORP effect. Loss of material from the sec¬ 
ondary stabilizes the remaining orbiting compo¬ 
nents. The lost mass may reaccrete onto the pri¬ 
mary, perhaps contributing to the observed equa¬ 
torial ridges, or may escape from the system. In 
these cases, the system undergoes another chaotic 
binary episode with three possible outcomes: a 
re-shaped asteroid, an asteroid pair, or a stable 
binary. These binaries still possess positive free 
energy such that they may disrupt if disturbed. In 
other cases, the system retains three components 
after secondary fission. While the numerical sim¬ 
ulations of Jacobson and Scheeres (2011b) did 
not yield this latter outcome, it is possible that 
this pathway explains the existence of stable triple 
systems. 

After stabilization of the low mass ratio bi¬ 
nary system, the secondary synchronizes due to 
tides (e.g., Goldreich and Sari 2009), although 
some satellites may be trapped in a chaotic ro¬ 
tation state for durations that exceed the classic 
spin synchronization timescales (Naidu and Mar¬ 
got 2015). Then the system evolves according to 
the BYORP effect and tides. These binary evo¬ 
lutionary processes and their outcomes are dis¬ 
cussed in Walsh & Jacobson (this volume). As 
shown in Fig. 11, these evolutionary paths include 
each of the binary morphologies identified in this 
chapter and by other teams (Pravec and Harris 
2007; Fang and Margot 2012c). In particular, the 
formation of wide asynchronous binaries such as 
(1509) Esclangona, (4674) Pauling, (17246) 2000 
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Rc-Shapcd Asteroid Asteroid Pair 


o 


Fig. 11.— Flowchart showing the possible evolutionary paths for an asteroid after it undergoes rotational 
fission. Each arrow is labeled with the dominant process and an estimated timescale for this process. Underlined 
states arc nominally stable for a YORP effect timescale. Figure from Jacobson and Scheeres (201 lb). 


GL 74 , and (22899) 1999 TO 14 is best explained by 
a rotational fission mechanism (Polishook et al. 
2011) followed by BYORP orbital expansion (Ja¬ 
cobson et al. 2014b). 

An alternative formation mechanism for triples 
such as (153591) 2001 SN 263 and (136617) 1994 
CC is that after creating a stable binary system, 
the primary undergoes rotational fission a second 
time. As long as the third component is on a dis¬ 
tant enough orbit, then this process may result in 
a stable triple system (Fang et al. 2011; Fang and 
Margot 2012c; Jacobson et al. 2014b). 


5. LARGE ASTEROIDS: SYNTHESIS 

The primaries of most known binary and triple 
asteroids greater than 20 km have spin periods in 
the range of 4 h to 7 h (Fig. 7). While these spin 
rates are not near the disruption spin limit, they 
are typically faster than the mean spin rates for as¬ 
teroids of similar sizes. The total angular momen¬ 
tum content, however, is well below that required 
for rotational fission. The secondary-to-primary 
mass ratios in these systems range from 10~ 6 to 
10 2 . These properties are consistent with satel- 
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lite formation during large collisions (Fig. 12). 
Durda et al. (2004) have shown in numerical sim¬ 
ulations that impacts of 10- to 30-km diameter 
projectiles striking at impact velocities between 
3 kms -1 and 7 kms” 1 can produce satellites that 
match observed properties. Multiple asteroid sys¬ 
tems, e.g., (45) Eugenia (Merline et al. 1999; 
Marchis et al. 2007) and (87) Sylvia (Margot and 
Brown 2001; Marchis et al. 2005a) can also plau¬ 
sibly form through collisions. 



Fig. 12.— Numerical simulations show that binaries 
can form as a result of large impacts between asteroids. 
In some scenarios, impact debris can remain gravita¬ 
tionally bound to the target body, forming a satellite 
(SMATs). This process likely explains the formation 
of large MBA binaries. In other scenarios, two frag¬ 
ments from the escaping ejecta have sufficiently sim¬ 
ilar trajectories, such that they become bound to one 
another (EEBs). Figure from Durda et al. (2004). 

There is more uncertainty related to the forma¬ 
tion of (90) Antiope and (617) Patroclus, which 
are both too large to be substantially affected by 
YORP. Hypotheses for the formation of (90) An¬ 
tiope include primordial fission due to excessive 
angular momentum (Pravec and Harris 2007), an 
improbable low-velocity collision of a large im- 
pactor (Weidenschilling et al. 2001), or shrinking 
of an initially wide binary formed by gravitational 
collapse (Nesvomy et al. 2010). Gravitational 
collapse in a gas-rich protoplanetary disk has 
been invoked to explain the formation of numer¬ 
ous binaries in the trans-Neptunian region. (617) 
Patroclus may be a primordial TNO that avoided 
disruption during emplacement in the trojan re¬ 


gion (Nesvorny et al. 2010). Wide TNO binaries 
would not be expected to survive this process, 
whereas encounter calculations (e.g., Fang and 
Margot 2012a) show that tight binaries would. 

6. CONCLUSIONS 

Studies of binaries, triples, and pairs remain 
a fertile ground for observing processes that are 
important in planet formation and for measur¬ 
ing quantities that are difficult to obtain by other 
means. These include masses and densities as 
well as thermal, mechanical, and interior proper¬ 
ties. Binaries or triples have been found in ~50 
NEAs, ~50 small MBAs, ~20 large MBAs, and 
2 trojans. A unifying paradigm based on rota¬ 
tional fission and post-fission dynamics explains 
the formation of small binaries, triples, and pairs. 
Because the sun-powered rotational fission pro¬ 
cess is unrelenting, and because the production 
of pairs is a frequent outcome of this process, a 
substantial fraction of small bodies likely origi¬ 
nated in a rotational disruption event. This origin 
affects the size distribution of asteroids and may 
explain the presence of single NEAs with equa¬ 
torial bulges observed with radar. Small satel¬ 
lites of large MBAs are likely formed during large 
collisions. Advances in instrumentation, observa¬ 
tional programs, and analysis techniques hold the 
promise of exciting findings in the next decade. 
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